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ABSTRACT 

We report the novel detection of complex high-column density tails in the probabil¬ 
ity distribution functions (PDFs) for three high-mass star-forming regions (CepOB3, 
MonR2, NGC6334), obtained from dust emission observed with Herschel. The low 
column density range can be fit with a lognormal distribution. A first power-law tail 
starts above an extinction (Ay) of ~6-14. It has a slope of a= 1.3-2 for the p oc r~ a 
profile for an equivalent density distribution (spherical or cylindrical geometry), and 
is thus consistent with free-fall gravitational collapse. Above Ay~40, 60, and 140, we 
detect an excess that can be fitted by a flatter power law tail with a >2. It correlates 
with the central regions of the cloud (ridges/hubs) of size ~1 pc and densities above 
10 4 cm -3 . This excess may be caused by physical processes that slow down collapse and 
reduce the flow of mass towards higher densities. Possible are: 1. rotation, which intro¬ 
duces an angular momentum barrier, 2. increasing optical depth and weaker cooling, 3. 
magnetic fields, 4. geometrical effects, and 5. protostellar feedback. The excess/second 
power-law tail is closely linked to high-mass star-formation though it does not imply 
a universal column density threshold for the formation of (high-mass) stars. 
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1 INTRODUCTION 

Probability distribution functions (PDFs) of column den¬ 
sity A(HI+H2) obtained using far-infrared emission of dust 


* E-mail: nicola.schneider@obs.u-bordeauxl.fr 


from Herschefi show a characteristic shape for low-mass 
star-forming regions: a lognormal distribution for low A, 
commonly attributed to turbulence, and a single power-law 

1 Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA. 
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tail for higher N (e.g., Schneider et al. 2013, and references 
therein). This sort of PDF is also found for extinction maps 
(Kainulainen et al. 2009; Froebrich & Rowles 2010). There 
is an ongoing discussion whether this power-law tail can be 
attributed to self-gravity (Klessen 2000; Kritsuk et al. 2011; 
Federrath & Klessen 2013; Girichidis et al. 2014; Schneider 
et al. 2013, 2015a) or is pressure-driven (Kainulainen et al. 
2011 ). 

In this Letter , we show that PDFs obtained from Her- 
schel dust column density maps of giant molecular clouds 
(GMCs) not only show a clear lognormal plus power-law 
tail distribution but can also exhibit an excess at very high 
column densities that can be identified as a second , shallower 
power law. It seems that this discovery is restricted to the 
densest regions of massive molecular clouds (ridges/hubs), 
though not all those GMCs display this feature. The objec¬ 
tive of this study is to report our detection and give some 
tentative explanations for its existance. 

2 OBSERVATIONS 

We make use of Herschel- derived column density ( N ) PDFs 
of three well-known GMCs. Data and PDFs of MonR2 at a 
distance of 0.83 kpc (Rayner et al., in prep, Didelon et al. 
2015) and NGC6334 at a distance of 1.35 kpc (Russeil et al. 
2013) were obtained within the HOBYS ( Herschel imaging 
survey of OB Young Stellar objects) key program (Motte 
et al. 2010), while the data for CepOB3 (distance 0.7 kpc, 
PI Gutermuth) were taken from the Herschel (Pilbratt et al. 
2010) archive. We selected these sources as examples because 
they cover different ranges in distance, mass, and UV-field 
(Table 1). Dust optical depths and temperatures were deter¬ 
mined from a greybody fit to the surface brightness at the 
160 pm wavelength of PACS (Poglitsch et al. 2010) together 
with the 250, 350, and 500 pm wavelengths of SPIRE (Grif¬ 
fin et al. 2010). For that, we assume a constant line-of-sight 
temperature for each pixel. The optical depths are then con¬ 
verted into H 2 column densities taking a dust opacity k\ = 
0.1 X (A/300pm)~P cm 2 /g with a fixed J3=2. Opacity varia¬ 
tions were observed in Orion A (Roy et al. 2014) so that we 
arrive to a relative accuracy of -50% over the whole column 
density range covered in Herschel observations. For more de¬ 
tails, we refer to Hill et al. (2011), Schneider et al. (2015a), 
and Konyves et al. (2015). For CepOB3, the PACS 160 pm 
data contain nonfunctional bolometers. The resulting col¬ 
umn density maps thus show a regular pattern with missing 
points, so that we fitted only the 250, 350, and 500 pm data 
points. The angular resolution of the maps is 36", matched 
to the longest-wavelength Herschel band. 

3 PROBABLITY DISTRIBUTION 

FUNCTIONS OF GMCS 

Figure 1 shows the column density maps and corresponding 
PDFs for the GMCs of our study. To describe the PDF, we 
use the notation p(^) with 

p(rj)dT] = (2 7tct^)~ 0 - 5 exp[-(;7 - fi) 2 / (2<t^)] dr\ (1) 

with Tj = In (A/(A)), cr^ as the dimensionless dispersion of the 
logarithmic field, and p as the mean. We determine cr r] with 
a fit to the assumed lognormal low column density part of 
the PDF and the slopes from a linear regression power- 
law fit with p(T]) = Po(t]/t]q) s to the high-density parts (see 


Schneider et al. 2015a,b, for more details). The intersec¬ 
tion between lognormal and power-law tail regime lies at 
extinctions 2 Ay of around 6-8 for CepOB3 and MonR2 and 
Ay—14 for NGC6334. Starting at these values, we then fit a 
power-law to the distribution up to the break point where 
additional high Ay excess becomes apparent, i.e., at Ay-40 
for MonR2, Ay-60 for CepOB3, and Ay-140 for NGC6334. 
This spread in values shows that the shift to a shallower 
power law does not represent a universal threshold such as 
the one proposed for massive star formation around Ay= 
300 (Krumholz & McKee 2008). The excess is clearly evi¬ 
dent over a range of Ay -40-200 for MonR2, Ay -60-200 
for CepOB3, and Ay -140-300 for NGC6334. Although it 
can be fitted by a power-law 3 , we cannot exclude the pos¬ 
sibility that other distributions can fit the excess similarly 
well. In considering its physical interpretation in comparison 
to the first power-law tail, we keep the premise of a second 
power-law. 

A word of caution is appropriate with regard to the PDF 
shape. Additional extinction due to foreground/background 
emission provokes that the width of the lognormal part be¬ 
comes more narrow and the peak shifts to higher Ay, while 
the slope(s) of the power-law tail(s) are steeper, compared to 
an undisturbed distribution (Schneider et al. 2015a). How¬ 
ever, the slope of the second power-law tail is less affected 
because a contamination of a few Ay does not significantly 
change the highest column densities. Contamination and 
projection effects, however, could still lead to substructure 
in the PDF, which might complicate the simple lognormal 
+ power-law tail(s) representation. The accuracy of the log¬ 
normal fit is thus limited and the first deviation point is not 
well defined. 

4 SPATIAL IDENTIFICATION OF THE 

DIFFERENT PDF REGIMES 

To identify the pixels that constitute the two power-law tails, 
we outline in Fig. 1 with a red and yellow contour the re¬ 
spective Ay-thresholds. While the red contour traces mainly 
filaments, the yellow contour encloses individual clumps with 
a size scale of -1 pc, delimiting single (MonR2) or several 
(CepOB3, NGC6334) star-forming center(s) of the GMCs. 
These central regions correspond in their geometry and 
physical properties to what was defined as ’ridges’ (or hubs) 
in massive clouds (Hill et al. 2011; Hennemann et al. 2012), 
i.e., regions with high column density (N >10 23 cm -2 ) over 
small areas. The PACS 70 pm images (Fig. 1) indicate the 
stellar/protostellar activity in these clumps: NGC6334 con¬ 
tains several Hu regions and a few ultra-compact (UC) Hu 
regions (e.g., Russeil et al. 2013), MonR2 hosts several UC 
H ii regions and bipolar outflows (e.g., Fuente et al. 2010; Di- 
erickx et al. 2015), and CepOB3 includes an embedded mas¬ 
sive cluster. Table 1 lists the masses contained above the two 
break points. For the high column density regime (yellow), 

2 For better comparison to other studies in the literature, we 
use the visual extinction value derived from the column density 
adopting N(H2)/Ay=0.94xl0 21 cm -2 mag -1 (Bohlin et al. 1978). 

3 The reduced x 2 rec i values of the second power law fits scat¬ 
ter around unity (1.05+0.28, 1.08+0.10, 2.43+0.16 for CepOB3, 
MonR2, NGC6334) indicating that the general assumption of a 
power-law dependence is justified. 
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Table 1. Properties of GMCs: total mass M derived above the first closed contour, i.e. Ay =2 for CepOB3 and MonR2 and Av=4 for 
NGC6334, masses (the percentage of the total mass given in parenthesis) above the first (DPI) and second (DP2) deviation point in 
the PDF (values determined from the Herschel N-maps), average and maximum UV-flux in units of the Habing field (obtained from the 
Herschel fluxes (Roccatagliata et al. 2013)), slopes sl/s2 of power-law tails, exponents 0 * 1/2 and o c i /2 for a spherical and cylindrical density 
distribution, equivalent beam deconvolved radius r (taking an area A=nr 2 ), and density n, determined from the average column density 
and r with n - N/(2r). The values of r and n for NGC6334 are averages over several clumps. The errors of both a c and a s determined from 
the first slope are -0.02, and those from the second slope are -0.05. 
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we determine also the equivalent radius r, and the density n. 
It becomes obvious that most of the mass is constituted by 
the lower column density regime, i.e., 83%, 61%, and 68% 
for CepOB3, MonR2, and NGC6334, respectively. The mass 
above the first deviation point makes up 15% to 39% of the 
complexes. The highest column density clumps still account 
for a significant amout of mass in MonR2 and NGC6334 
(both 15%), but a much smaller proportion (<1%) for Ce- 
pOB3. The average density (n) in these clumps is also high, 
i.e., between 2xl0 4 cm -3 and 3.5xl0 5 cm -3 . However, the 
existence of high density clumps/cores does not by itself ex¬ 
plain the second power law. In a region such as Aquila, ~15% 
of the total cloud mass resides in gas of density larger than 
10 4 cm -3 (Konyves et al. 2015), but only a single power-law 
tail is observed (Schneider et al. 2013; Andre et al. 2014). 
We focus now on the analysis of the two power-law tails. 
From the slope s, the exponent a for the density distri¬ 
bution p(r) oc r~ a can be derived. For a spherical geome¬ 
try, representing clumps and cores, a = -2/s+ 1 (Federrath 
V Klessen 2013). For a cylindrical one, characterizing fila¬ 
ments, a = -I/ 5 + 1 (Myers 2015). Table 1 lists the values for 
both geometries for the first and second tail. The values of 
a for the first tail range between 1.3 and 2, consistent with 
free-fall collapse, see Schneider et al. (2013, 2015a, b) and 
Girichidis et al. (2014) for a detailed discussion of both fila¬ 
ments and clumps/cores. The values of a c 2 between 1.85 and 
2.64 correspond very well to the range of values directly mea¬ 
sured for the B211/B213 filament in Taurus in Palmeirim et 
al. (2013). Independent observational support for gravita¬ 
tional contraction comes from molecular line observations 
of other GMCs that show spectral infall signatures across 
ridges (Schneider et al. 2010, 2015b; Galvan-Madrid et al. 
2010; Peretto et al. 2013). These studies show that the for¬ 
mation of high-mass stars requires a very large mass ac¬ 
cretion rate, provided by infall from merging filaments and 
gravitational collapse of larger structures on parsec-scales. 
For the excess - presuming it can be described by a power- 
law - a lies between 1.9 and 4.2. The assumption of a spheri¬ 
cal density distribution is probably more valid here, because 
the corresponding regions are more compact and circular. 
However, free-fall alone can not produce a >2, only a slowed- 
down collapse can lead to piling up high (column) densities 
on very small spatial scales. 


5 SIGNIFICANCE OF THE EXCESS 

The PDF excess we observe for the GMCs in our study is 
not a systematic property of all GMCs. It is also seen in 
Rosette (Schneider et al. 2012), W43 and W51 (Schneider 
et al., in prep.), and W3 (Rivera-Ingraham et al., 2015), but 
there is no example of a low-mass star-forming cloud that 
displays this feature. Similarly, there are also massive clouds 
such as Carina and NGC3603 (Schneider et al. 2015a), Vela 
(Hill et al. 2011), or M16 (Hill et al. 2012) where no clear 
indication for this excess is found in the PDFs of the whole 
cloud. Therefore, it is unlikely that the excess is caused by a 
systematic bias introduced during the determination of the 
dust column density maps. 

It is out of the scope of this paper to discuss in detail 
the general uncertainties related to the determination of the 
column density maps and we refer the reader to Roy et al. 
(2013, 2014) and Konyves et al. (2015). Nevertheless, we 
mention here briefly the biases that could potentially influ¬ 
ence the inferred slope and are thus relevant for the results of 
this paper: 1) An overestimation of the line-of-sight temper¬ 
ature (assumed to be constant for the SED fit) leads to an 
underestimation of the column density (and the vice versa). 
2) If dust opacity increases towards high column densities, 
then this could give rise to a steeper slope of the power- 
law tail at higher densities. 3) The adopted value of the 
specific dust opacity /3—2 could be too high in parts of the 
warm GMCs. A lower (higher) value of ft in the greybody fits 
leads to higher temperatures and thus lower (higher) column 
densities. 

6 WHAT CAUSES THE TWO POWER-LAW 

TAILS IN THE PDF? 

Given that the excess/second power-law tail is only found in 
high-density (( n ) ~10 4-5 cm -3 ) star-forming clumps on size 
scales of ~1 pc, i.e. in ridges/hubs, the physical process(es) 
in control of it must be active on small scales. It is evident 
that at these densities gravitational collapse of individual 
star-forming cores is involved, but it is not yet understood 
what produces the remarkable excess in column density. 

We emphasize that in all cases the newly found excess 
is preceeded by a first power-law in the PDF consistent with 
the dominant effect of self-gravity found by numerical sim¬ 
ulations (e.g., Vazquez-Semadeni et al. 2008; Kritsuk et al. 
2011) as well as observations (e.g., Schneider et al. 2013, 
2015a,b). We thus expect the regions under consideration to 
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Figure 1. Left: Dust column density maps of the GMCs. Contour lines in red and yellow indicate the break points between the 
lognormal distribution and the first and second power-law tail, respectively. The small panels show a zoom into the central regions, 
displaying PACS 70 /im emission (scale typically 0.1-30 Jy/pixel). Right: Corresponding PDFs at an angular resolution of 36" (binsize 
0.1). The left y-axis indicates the normalized probability p{rf} : the right y-axis the number of pixels per logarithmic bin. The error-bars 
were calculated using Poisson-statistics. The upper x-axis gives the extinction Ay and the lower x-axis 77 = ln(N/(N)). The green curve 
indicates the fitted lognormal PDF and the red and blue lines the power law fits to the high-density tails. The fit to the first tail was 
only performed between the break points but the line was continued to emphasize the excess. The deviation points (DP) and slopes (s 
with its error) of the power-law tails are given in each panel. 

contract in almost free-fall (Girichidis et al. 2014), i.e., gas 
at a density pi falls towards higher densities p' > p;. Without 
any change in physical processes for higher densities, free-fall 
would proceed unimpeded and the power-law would extend 
indefinitely. Deviations from this idealised single power-law 
can be twofold. On the one hand the second excess could 


be due to a change in dynamics at high densities or due to 
observational effects. In the following we give some tentative 
explanations that, however, need more profound studies. 

Any physical process that slows down the free-fall mo¬ 
tions reduces the flow of mass towards higher densities. We 
stress that in this paradigm the deceleration of free-fall stems 
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from within the centres of the cloud clumps. Processes act¬ 
ing from the outside like a hot ambient medium in which 
the cloud is embedded or large-scale colliding flows are un¬ 
likely to directly alter the PDF at the highest densities with¬ 
out altering the distribution at lower densities that is in the 
free-fall regime of the first power-law. Rotational effects (1.) 
have been invoked to explain excess at high densities seen in 
a density PDF obtained from an isothermal, self-gravitating 
supersonic turbulence simulation (Kritsuk et al. 2011). The 
spatial scales at which angular momentum is likely to dom¬ 
inate, however, are probably smaller than the resolution of 
the observations presented in this paper (0.1-0.25 pc). Like¬ 
wise, thermodynamical effects of increasing thermal pressure 
due to shielding and reduced cooling (2.) is unlikely to dom¬ 
inate at the observed scales (e.g., Larson 2005). The role 
of magnetic fields (3.) is rather unclear. It was shown that 
strong (Koertgen et al. 2015) or intermediate (Heitsch et al. 
2001 ) magnetic fields, acting on a clump scale, slow down 
and can even completely prevent the star-formation process 
in magneto-hydrodynamic (MHD) simulations. If magnetic 
fields entirely prevent the collapse and star formation, the 
clouds are stable and will not show a first gravitationally 
dominated power-law in the PDF. Moderate and weak mag¬ 
netic fields tend to reduce the degree of fragmentation but 
once a gravitational instability sets in, the fields are unlikely 
to stop further collapse and star formation (e.g., Ziegler 
2005; Banerjee et al. 2006). In this case, the N-PDF shows 
a power-law tail due to the gravitational collapse of the su¬ 
percritical cores. Moreover, a change in dominant geometry 
(4.), i.e., from filamentary to spherical, could also provide an 
explanation. Longitudinal filament collapse on parsec scales 
reduces the mass transfer rate and dense gas is then piled 
up (Toala et al. 2012). 

All of the above proposed mechanisms to produce col¬ 
umn density excess should apply for all types of clouds - 
not only to the most massive GMCs. Therefore, high col¬ 
umn density excesses should be observed more commonly, 
but this is not the case. Note, however, that the slope of the 
power-law tail(s) depends on the projection, i.e. the view¬ 
ing angle, in which the cloud is being observed (Ballesteros- 
Paredes et al. 2011). In any case, a major difference be¬ 
tween clouds forming only low-mass stars and those with 
high-mass star-formation is stellar feedback (5.). For exam¬ 
ple, additional compression by internal ionization due to an 
ultra-compact H n region can provoke a power-law tail with 
a >2 (Tremblin et al. 2014). (High-mass) outflows may also 
play a role. It was shown by Sadavoy et al. (2014) that the 
slope of the power-law PDF in different regions in Perseus 
depends on the local feedback from low-mass young stellar 
objects. A recent study of MonR2 (Dierickx et al. 2015) re¬ 
vealed a significant number of CO outflows in the central 
region that imprints on the velocity structure of the region. 
Not all clouds with protostellar feedback, however, also show 
a second power-law. If feedback does not efficiently couple 
to the surroundings and the dense gas, an excess in the PDF 
might not develop or might not be visible at the observed 
scales in the column density PDF. 

7 CONCLUSIONS AND OUTLOOK 

At this point, it is not possible to give a final answer for our 
detection of excess/second power-law in the PDF of mas¬ 
sive GMCs. Further studies are required to look into more 


detail which of the proposed processes can play a dominant 
role. Besides more sophisticated 3D hydrodynamic simual- 
tions of star-forming regions a deeper understanding of how 
local thermal and dynamical properties are reflected in ob¬ 
servations is needed. 
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